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We have searched for prompt production of Xci, Xc2 and X(3872) in continuum e^e~ annihilations 
using a 386 fb~^ data sample collected around y/s — 10.6 GeV with the BABAR detector using the 
'yJ/'4> decay mode. After accounting for the feed-down from ^l'{2S) 7Xci,2, no significant signal 
for prompt Xci,2 production is observed. We present improved upper limits on the cross-section, 
with the rest of the event consisting of more than two charged tracks, to be 77 fb for Xci and 79 fb 
for Xc2 with e"'"e~ center-of-mass frame Xc momentum greater than 2.0 GeV at 90% confidence level. 
These limits are consistent with NRQCD predictions. We also set an upper limit on the prompt 
production of X(3872) through the decay X(3872) ^ jJ/ip. 

PACS numbers: 13.66.Bc, 12.38.Qk, 12.38.Bx, 14.40.Gx 



Charmonium production in e~^e~ annihilation pro- 
vides opportunities to study both perturbative and non- 
perturbative effects in QCD and to search for new char- 
monium states [l|, 0]. The prompt production of J/^p 
and ■0(25') in e"'"e~ annihilation i^] and of double- 
charmonium 0, 0] have been observed at i?-factory ex- 
periments. These observations are surprising because the 
measured cross-sections are larger than non-relativistic 
QCD (NRQCD) calculations by up to an order of mag- 
nitude [l|,|7l. 

In the NRQCD production mechanism, a heavy 
quarkonium (qq) state can be produced at short dis- 
tances as a conventional color-singlet, or as a color-octet 
state, which then evolves into an observed quarkonium 
meson along with other light hadrons. With this color- 
octet mechanism, one may explain the enhancement for 
J/ijj production in e+e~ annihilation [ij. The produc- 
tion of Xci,2 (Xc) in e+e^ annihilations is an excellent 
probe of color-octet contributions, which are more promi- 
nent in Xc production than in J/tJj production. This is 
because color-octet and color-singlet processes enter Xc 
production at the same order, and C-parity suppresses 
the process e+e^ ccgg, which dominates J/ip produc- 
tion. Calculated cross-sections for prompt Xc produc- 
tion in e^e^ annihilations are a{e^e^ — > Xci^) = 85 
fb and cr(e"*"e~ — > Xc^X) = 123 fb, with e^e^ center- 
of-mass (CM) frame J/^ momentum p}/^ > 2.0 GeV, 
and where X is one of qq, gg and g in the leading order 
processes and J/V' is from Xc decay Q. More 

accurate measurements for Xc states will help to clarify 
the discrepancy between theoretical calculations and ex- 
isting measurements, and may point to other methods 
and mechanisms in QCD to explain the differences. 

Prompt production of charmonium mesons in e+e~ an- 
nihilation has been searched for using either the recon- 
structed mass in an exclusive decay mode 0, 01 or the 
mass distribution of the system recoiling against the J/ip 
or ^{2S) 1,0. Although prompt production of XcO has 
been observed, prompt production of the other Xc states, 
Xci and Xc2, has not been observed. 

In this paper, we present a search for prompt Xc pro- 
duction in continuum e~^e~ annihilation using the jJ/ip 
{J — > decay mode, which is experimentally clean 

and is the dominant one in Xc decay. The current limits 
on prompt production of Xc are a{e^e'' — s- Xci-''^) < 350 



fb and cr(e+e- Xc2^) < 660 fb with e+e" CM frame 
Xc momentum p* ^ > 2.0 GeV, where X is the rest of the 
event [1] . Belle and BABAR recently observed an indica- 
tion of the decay X(3872) li> in B decays [o'l, and 
therefore we also search for prompt X(3872) production 
using the ^J/^ decay mode in e^e~ annihilation. 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric-energy e+e" 
collider, where 9.0 GeV electrons and 3.1 GeV positrons 
are collided at a CM energy of 10.58 GeV, the mass of the 
T(4S') resonance. The integrated luminosity (£) consists 
of 349 fb-i {Con) at the T(4S') resonance and 37 fb'^ 
{Coff) at a center-of-mass energy 40 MeV below the res- 
onance. 

The BABAR detector is described elsewhere [l3| and 
here we give only a brief overview. The momenta of 
charged particles are measured by the silicon vertex 
tracker, consisting of five layers of double-sided silicon 
strip sensors, and the central drift chamber (DCH) with 
40 wire layers, both operating in a 1.5 T magnetic field 
of a solenoid. The tracking system covers 92% of the 
solid angle in the CM frame. An internally-reflecting 
ring-imaging Cherenkov detector (DIRC) with quartz bar 
radiators provides charged particle identification (PID). 
A CsI(Tl) electromagnetic calorimeter (EMC) is used to 
detect and identify photons and electrons, while muons 
are identified in the instrumented magnetic flux return 
system (IFR). 

Electron candidates are identified by the ratio of the 
shower energy measured in the EMC to the track momen- 
tum measured in the DCH, the shower shape, the specific 
ionization energy loss in the DCH, and the Cherenkov an- 
gle measured by the DIRC. Muons are identified by the 
depth of penetration into the IFR, the IFR cluster ge- 
ometry, and the energy deposited in the EMC. Photon 
candidates are identified by EMC clusters that have a 
shape consistent with an electromagnetic shower and are 
not associated with a charged track. 

We use a Monte Carlo (MC) simulation of the BABAR 
detector based on GEANT4 [llj to validate the anal- 
ysis procedure, to evaluate signal detection efficiencies, 
to model probability density functions (PDFs), and to 
estimate background contributions. We use samples of 
e+e~ Xci, 2 + J/'4' OY ■0(25') MC events to deter- 
mine the selection criteria. To estimate the signal re- 
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construction efRciencies and PDFs, we use single Xc 
MC samples decaying to ^Jj'ii with J/V" e+e" or 
Jj^l) , which are generated with flat distribu- 

tions in p* (CM frame Xc momentum) and cos 9* (cosine 
of the polar angle of the Xc momentum to the beam axis 
in the CM frame). To understand combinatorial back- 
ground, we use MC generated e"'"e~ — *■ rjcXcO, or ric{2S) 
events produced in association with either J/ip or ipi^S) 
mesons. BB generic and initial state radiation (ISR) 
ip{2S) {e~^e~ jip{2S)) MC events are used to estimate 
background contamination. The Xc candidates from B 
decay are used as a control data sample to correct for 
differences in the photon energy measurements between 
MC simulation and data. 

Charged particles are required to have a point of clos- 
est approach to the beam spot of less than 10 cm along 
the beam axis and less than 1.5 cm in the plane trans- 
verse to the beam. The J/tp mesons are reconstructed in 
the dilepton channel using two oppositely-charged tracks 
identified as electrons or muons. An algorithm to recover 
the energy loss due to bremsstrahlung is applied to elec- 
tron candidates. The invariant mass of the reconstructed 
J/ip is required to be within the range [3.07, 3.13] GeV 
for the fi^n^ channel and [3.05, 3.13] GeV for the e+e^ 
channel. The asymmetric selection in the e^e^ chan- 
nel is due to initial and final state radiation. The J/tp 
candidate is subjected to a vertex constrained fit and is 
combined with a photon candidate that satisfies stan- 
dard reconstruction quality criteria as described below. 
Multiple signal candidates in the event are allowed. 

The photon candidates are EMC clusters in the angu- 
lar region 0.41 < 9 < 2.41 radians where 9 is the polar 
angle with respect to the beam axis in the laboratory 
frame. The lateral energy distribution {LAT) [1^ mea- 
sures the transverse energy profile of a cluster; requir- 
ing this to be less than 0.5 suppresses clusters due to 
both electronic noise and hadronic interactions. The az- 
imuthal asymmetry of the energy deposition in a cluster 
is measured by the A42 Zernike moment [l^. Requiring 
A42 less than 0.1 further rejects clusters from hadronic 
interactions. In addition, the angular separation between 
the direction of the candidate and of any charged track 
in the event should be at least 9° in the laboratory frame 
(spht-off rejection). The clusters satisfying these criteria 
come mostly from tt" decay. We reject photon candidates 
that, when combined with any other photon, produce a 
mass between 114 MeV and 146 MeV (7r° veto). The 
partner photon must have energy greater than 30 MeV 
and LAT < 0.8 without any requirement on A42 and 
spht-off rejection. 

Backgrounds arise from combinatorial background in 
B decays and continuum events, and decays of ipi^S) 
mesons produced either promptly or in ISR events. To 
suppress B-background contributions, we require p*^ > 
2.0 GeV and > 2.0 GeV. For the Xc control sam- 

ple, we require p*^ < 1.7 GeV and < 2.0 GeV. The 



combinatorial background for J/tp candidates in contin- 
uum events is reduced by requiring \cos9'^^\ < 0.9 0], 
where 9'^^ is the J/ip helicity angle, measured in the rest 
frame of the J/ip, between the positively charged lepton 
daughter and the jJ/tp system. 

The backgrounds from prompt ipi^S) radiative decay 
to are indistinguishable from the signal. The esti- 
mated contribution from prompt ipi^S) production will 
be subtracted from the measured cross-sections. 

Substantial backgrounds are due to ISR production of 
'0(25') decaying to 7Xc, which produces low multiplicity 
and a jet-like event shape. To suppress such backgrounds, 
the ratio R2 of second and zeroth Fox- Wolfram moments 
of the event [3] is required to be less than 0.8 and the 
number of charged particles in the event is required to 
be at least five (Nch cut). We estimate the possible con- 
tributions from ISR production using MC samples and 
subtract them from the signal yield. Two-photon back- 
ground contributions are estimated to be negligible with 
all selection criteria applied. 

The helicity angle of the jJ/ip system {9h) is the angle, 
measured in the rest frame of the ^J/ip system, between 
the momentum of the J /ip and the momentum of the 
e+e^ center-of-mass in the laboratory frame. The J /ip 
mesons from combinatorial background tend to be along 
the direction of the boost vector which makes cos 9h close 
to unity whereas the distribution of signal events is flat. 
We optimize the cos 9h cut using MC samples by max- 
imizing the figure of merit N'^^/ [Ncont + ^bb)' where 
Nsig, Ncoiit and Ng-g are the numbers of events from 
signal, continuum, and BB background expected in the 
data sample respectively. The scale of Nsig is not sen- 
sitive to the optimized cut. For Ncont we use the yield 
from off-resonance data multiplied by (Con+^of f)/ ^off- 
The optimized cut is found to be cos 9h < 0.4. The same 
cut is applied for the X(3872) search which has similar 
kinematics. 

We extract the signal yield using an unbinned maxi- 
mum likelihood (UML) fit (nominal fit) for the distribu- 
tion of AM, the mass difference between the signal Xc 
or X{3872) candidate and the daughter J/tp candidate. 
We use a AM range [0.25, 0.60] GeV for the Xc searches 
and [0.60, 0.95] GeV for the X{3872) search in the nom- 
inal fit. To estimate the systematic uncertainty, we use 
[0.25, 0.35] GeV and [0.50, 0.60] GeV as sideband regions 
and [0.35, 0.50] GeV as the core signal region for the Xc 
states. 

The AM distribution for signal candidates is described 
by a Crystal Ball Line shape (CBL) which is a Gaussian 
(described by the peak value AMq and resolution a am) 
with a power law tail l/( AMq — AM -I- const)", at a value 
of AMq — a ■ a AM- We use different PDFs for Xci, Xc2, 
and X(3872), averaged over the e+e" and ^^/^^ modes. 
The parameter values used in the CBL are determined 
using MC simulation and are then fixed in the nominal 
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To estimate the signal detection efficiency e, we de- 
compose it into three factors: efficiencies of reconstruc- 
tion (e^), tt" veto (£„) and spUt-off rejection {eg). The 
efficiency becomes smaller in low p* bins and high cos 9* 
bins owing to the P*j/^ > 2.0 GeV requirement and lower 
detector coverage near the endcap region. To get an esti- 
mate of Cr, we divide the region 2.0 < p* < 5.0 GeV into 
6 bins and —1.0 < cos 6** < 1.0 into 5 bins. We correct 
using the formula — wf Wj°^^ where wf and 
^cosS weights and e.y is an efficiency matrix {i —1,6; 
j =1,5), averaged over the e'^e~ and pL~ modes using 
single particle MC samples. The weights are defined by 




0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 
AM (GeV) 



FIG. 1; (color online) Nominal fit result for the Xc search 
(top) and the X(3872) search (bottom) from the 386 fb~^ 
data sample (p* > 2.0 GeV). The points represent the data, 
the dashed lines are background PDFs, the solid lines below 
the points are signal PDFs, and the solid lines on the points 
are the total PDFs. 



fit. The resolution ctam is 14.0 MeV, 15.3 MeV and 20.5 
MeV for Xci, Xc2, and X(3872) respectively and these 
are scaled by /?, a scale factor for the AM resolution. 
The mean AMq for each of xa, Xc2, and X(3872) is 
given by the known mass shifted by S, an offset of the 
PDF in AM. The difference of the Xci and Xc2 masses 
is constrained to the known value 45.5 MeV [l3|. The 
P and S parameters are determined as /3 = 0.89 ± 0.03 
and 6 = (2.7 ± 0.4) MeV using a control data sample of 
Xc mesons from B decay and fixed in the nominal fit. 
The background line shape is described by a third-order 
Chebyshev polynomial with free coefficients. 

The results for the nominal fit are presented in Fig. 
[TJ For the Xci and Xc2 searches, we analyze 1417 events 
after all selection criteria. The number of Xci candidates 
is 134 ± 23 and the number of Xc2 candidates is 56 ± 19. 
For the X(3872) search, we find iVx(3872) = -8± 11 from 
293 events. 

The ISR ip{2S) backgrounds are estimated using MC 
samples to be 9.4 events for Xci and 5.1 events for Xc2- 
Subtracting these from the fitted yields we find N^^-^ ~ 
125 ± 23 and N^^^ = 51 ± 19, which we attribute to 
the sum of prompt Xc production and feed-down from 
prompt ■0(25') production. 
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cos 0* 



where and e™"^ are efficiencies in bins of p* and 
cos^*, determined from the single Xc MC samples, and 
and are the yields in each bin, extracted 

from the binned fit to the data sample. For the X(3872) 
search, we use the averaged efficiency when the weights 
for Xci and Xc2 are used, because of the limited statistics 
for the number of X(3872) candidates. The Cij values are 
determined from the single X{3872) MC sample. With 
these corrections, the values are 10.1%, 9.3%, and 8.4% 
for Xci, Xc2, and X(3872), respectively. 

To estimate £„ and e^, we need to have knowledge of the 
efficiency as a function of photon (iV^) or charged track 
multiplicity (Nch), and the or Nch fractional distri- 
bution of signal events, because and es are strongly 
dependent on the number of photons or charged tracks 
in the event. We estimate efficiencies for each A^^ and Nch 
bin using signal MC simulation corrected by the data-to- 
MC difference using Xc candidates from B decays. The 
distributions of Nj and Nch for signal events are esti- 
mated from the sideband-subtracted data sample. The 
Nj distribution ranges from 1 to 18 and the Nch distri- 
bution ranges from 5 to 14. We estimate e^— 0.80 and 
es= 0.96 from an average calculated by the following for- 
mula: 



where Npi stands for the number of photons or charged 
tracks produced in the ith bin, No 



photons or charged tracks observed in the ith bin, ivpj 



for the number of 

Np^ = 

Noi/e{Ni), and e{Ni) is the efficiency of the ith bin in the 
distribution of N^ or Nch- For the X(3872) search, we 
use the same e„ and as for the Xc- The total efficiency e 
is the product Er- Cs and is estimated to be 7.7%, 7.1% 
and 6.4%, respectively, for the Xci, Xc2, and X{3872). 
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TABLE I: Systematic uncertainties (quoted in %) on (tjv^^>2 
defined in the text. 







Xc2 


A (3872) 


p*/cos6* correction 


13.3 


26.5 


34.9 


iracit euiciency 


U.O 


U.O 


U.O 




7.2 


7.2 


7.2 


Photon FID 


1.8 


1.8 


1.8 


tt" veto efficiency (ei,) 


2.3 


2.3 


2.3 


Split-Off rejection efficiency (ej,) 


0.4 


0.4 


0.4 


PDF 


3.5 


11.2 


15.1 


ISR Baclcground 


3.8 


5.0 






5.4 


5.0 


0.7 


Total 


17.1 


30.6 


38.8 



The sources of systematic uncertainty are summarized 
in Tabic [D The dominant uncertainty is from the recon- 
struction efficiency (e^) correction from the p* and cos 6* 
distributions. For the Xc search, we assign the system- 
atic uncertainty as the r.m.s. spread of 10,000 simulated 
experiments (each experiment gives one value) with 
weights generated according to the central values and er- 
rors from the p* and cos6'* binned fit results. For the 
X(3872) search, we adopt a conservative approach. We 
calculate separately the r.m.s. values corresponding to 
the binned fit results for Xci and Xc2, and assign the sum 
of r.m.s. values as the systematic uncertainty for the 
X(3872) reconstruction efficiency. 

The error from the PDF modeling is estimated by a 
quadratic sum over the changes in the yield from an al- 
ternative background line shape e-[Po+Pi('^*^)+P2(^^^)'] , 
and ±1 standard deviation of the uncertainties in the 
measured 5 and /3 in the Xc control sample from B decay. 
We take the data-to-MC difference in track reconstruc- 
tion efficiency as a source of systematic uncertainty. To 
estimate systematic uncertainties in charged PID efficien- 
cies, we assign the difference when taking ±1 standard 
deviation of each error depending on momentum and az- 
imuthal angle of tracks measured using control samples. 
The systematic uncertainty of photon identification is es- 
timated by comparing data with MC simulations of r+ — s- 
■K^v and r+ p{-k^tt^)v samples. We assign half of 
the ISR ip{2S) background estimate as systematic uncer- 
tainty for the Xc search. The uncertainty of the ISR back- 
ground is neglected for the X(3872) search. The JJBi is 
a product of sub-decay mode branching fractions, that is 
BiXc ^ iJ/^) ■ m-J/^P ^ e+e-) + S(J/V' ^ /i+A*-)]. 
The systematic error related to JJBi is estimated from 
the reference values [11]. The systematic uncertainties 
from e„ and Es evaluations are estimated by a quadratic 
sum over the deviations in two cases: when the data-to- 
MC correction is not used and when and Nch dis- 
tributions are taken without sideband subtraction to see 
the effect of backgrounds on the distribution. 

Table Hi] summarizes the measurements and all the 
quantities we need to calculate crjv^^>2, that is the cross- 
section of prompt Xc or X{3872) production {(7{e'^e~ 



ccX)) times the probability of the rest of the event {X) 
having more than two charged tracks, pN^h>^- The 
result crjVeh>2 is derived from the formula Ngg = £ ■ 



is the number of Xc or 
annihilation. In the case 



e ■ o-Ar<:h>2 • UBi where Nsg 
X{3872) candidates from e+e 
of Xc <7N^h>2 includes the prompt ip{2S) feed-down con- 
tribution. For the X(3872), we measure the product 
<JN^u>2 ■ B{X{i872) -fJ/xP) because the X(3872) ^ 
7J/V' BF is unknown. 

TABLE II: Signal yield Nsg from the nominal fit after sub- 
tracting the ISR tp{2S) estimate; signal detection efficiency 
(e = er ■ tv ■ Es); product of sub-decay mode BF's (H'^O; inte- 
grated on- and off-resonance luminosity {£,); crjv^^>2 (defined 
in the text) and its upper limit including systematic uncer- 
tainties; cr^°'>2' (''"^ch>2 f'^'' ^ii^ prompt production) and its 
upper limit including systematic uncertainties. Upper limits 
are at the 90% C.L. Note that ajv ^>2 for X(3872) denotes 
<Tiv,,>2-B(X(3872) ^7J/V'). 





Xcl 


Xc2 


X(3872) 




125±23 


51±19 


-8.0±11 


e{%) 
UB. (%) 
C (fb-i) 


7.7 
4.2 


« 75) 
7.1 
2.4 


(<15) 
6.4 
11.9 


386 


386 


386 


Cr]V,^>2 (fb) 


99±18 ± 17 


78±28 ± 24 
< 125 


-3±4 ± 1 
< 5 


<T>2 (fb) 


41±18±21 


23 ± 28 ± 26 


-3±4 ± 1 




< 77 


< 79 


< 5 



For prompt Xc production, it is necessary to sub- 
tract prompt ip(2S) feed-down to Xc- The contribu- 
tion of prompt ip{2S) production is estimated to be 
(58 ± 12) fb for Xcl and (54 ± 11) fb for Xc2 us- 
ing cr(e+e~ ^ ■^{2S)X) = (0.67 ± 0.13) pb for p* > 



151. The 



2.0 GeV 14] and the ip{2S) -fXc BF 
errors are included as systematic uncertainties in the 
prompt Xc production cross-section. Feed-down from 
other ip{2S) decay modes with photons is checked us- 
ing MC simulation: J / ipn'^ n'^ , J/ipr]{jj), J/iprjljir^iT'^), 
J/ipri{Tr'^Tr~TT^), and J/iljr]{TT°TT^n^). No background 
from these decays is seen in the MC simulation. The 
resultant cross-sections, (^^°^2 ^ ^'^^ production are 
shown in Table [Til 

Our measurements use an additional kinematic cut 
p*^ > 2.0 GeV which has little effect on the cross-section 
because leading-order contributions are from two-body 
e''~e~ annihilation processes. To compare these results 
with the theoretical predictions in Ref. the value 
of VNah>2 should be estimated correctly. Nevertheless, 
our upper limits are comparable with the NRQCD cross- 
section predictions. 

In summary, we have searched for prompt production 
of Xcl and Xc2 in e+e^ annihilation near ^/s = 10.6 
GeV. We observe candidates for these Xc states, but the 
measured cross-sections are compatible, within statistics. 
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with the expected contributions of Xc feed-down from 
prompt iIj{2S) production. The 90% confidence level 
upper limits on cr]^°™2* "^"^ ^ ^'^^ Xci and 79 fb 
for Xc2 with p*^ > 2.0 GeV. We find no evidence for 
prompt X(3872) production via the decay X(3872) 
jJ/tp- We set the 90% confidence lever upper limit on 
<7N,H>2 ■ B{X{3872) 7J/V') to be 5 fb. The upper lim- 
its presented on prompt Xc production are significant im- 
provements on the previously reported results [1| . These 
limits are comparable to the theoretical cross-section pre- 
dictions of Ref. [H . Upper limits on prompt production 
of Xc in comparison with J/ip and tpi^S) prompt produc- 
tion 0, Q , can be used to further our understanding of 
the charmonium prompt production mechanism [l|, U Q . 
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